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Two decades ago, Sadovskii found an exact solution of a model describing 
a pseudogap in electron energy spectrum (first introduced by Lee, Rice and 
Anderson). The discovery of a pseudogap in high-T c superconductors has 
revived the interest to his exact solution. I review the model with the emphasis 
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. on physical content, point out an error in the original Sadovskii's solution 

>> \ ! and explain which problem he actually solved. A recent incorporation of 

Sadovskii's ideas into a description of "hot spots" on the Fermi surface in 
cuprate superconductors (Schmalian, Pines and Stojkovic) is briefly discussed. 
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I. INTRODUCTION 



A model of electrons with a pseudogap from fluctuations of an order parameter was 
introduced in 1973 by P. A. Lee, T. M. Rice and P. W. Anderson [|IJ. A few years later, 
M. V. Sadovskii showed that it admits an exact solution The model describes a Peierls 
system (a metallic chain with a charge density wave instability) above the phase transition 
temperature Tp. The exact solvability comes at a price: (a) The solution is specifically 
tailored for one dimension, (b) It is assumed that Peierls-Kohn phonons are described by a 
non-selfinteracting boson field. These two limitations of the Sadovskii's solution have been 
known since its publication. 

Recently, however, I discovered an unfortunate error in the original paper by Sadovskii 
and now I am convinced that he actually solved a completely different, rather unphysical 
problem. This and the fact that Sadovskii's work is often regarded as the one and only exact 
model of the pseudogap |||| has prompted me to review this model. While its mathematical 
side has been discussed quite thoroughly by Sadovskii himself, the physical content deserves 
further comment. 

The plan of the paper is as follows. After a brief description of the Peierls instability in 
a one-dimensional conductor (Sec. |1 A|), a suitable mathematical formalism will be presented 



in Sec. fB]. It will be shown that assessment of multi-phonon contributions to the fermion 
energy spectrum requires a knowledge of statistical properties of the phonon ensemble. The 
model of Sadovskii, which postulates Gaussian statistics for the phonons, is introduced, 
interpreted and thoroughly illustrated in Sec. O. This is done in order to demystify its well- 
known yet strange-looking electron spectrum in the limit of long-range phonon correlation 
length £. I will then point out a previously unnoticed error in Sadovskii's "exact" solution 
for a finite £ (Sec. PTT|) and explain which problem Sadovskii has actually solved. Finally, a 
recent extension of the Sadovskii model to "higher" dimensions by Schmalian et al. || in 
the context of high-T c superconductivity will be discussed in Sec. [TV[ 



*Current address: School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540. 
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A. Peierls instability and fluctuations in 1 dimension 

An excellent introduction to the Peierls effect can be found in G. Griiner's book |J. See 
also an article by G. A. Toombs 0j, which reviews in detail theoretical and experimental 
developments prior to 1984. In a one-dimensional electron gas (with Fermi momentum pp) 
fluctuations of electron density are particularly strong near the wavevectors ±2pp. This 
happens because creating a hole and an electron with momenta near ±pp costs little energy. 
Therefore, when electron-phonon interaction couples lattice vibrations to these fluctuations, 
phonon modes with momenta near 2pp become "soft" and a static charge- density wave 
(CDW) appears below a transition temperature T P . 

In the context of the mean-field theory of a Peierls transition, an energy gap opens at 
the Fermi points exactly at T — Tp. However, some remnant of the gap can be created 
by fluctuations even above T P . Lee, Rice and Anderson calculated the electron self-energy 
induced by the emission and reabsorption of a (dressed) phonon (Fig. [[]). Approximating 
the two-point phonon correlation by two Lorentzian lines peaked at ±2pp, one obtains |1[ 

5 2 

Z(uj,pp + p) = — ■ v (1) 

to + pv + ivt, 

where S 2 = 47i 3 T P ^/7((3)v can be regarded as the average fluctuation of the order parameter, 
£ is the phonon correlation length and v the Fermi velocity. 
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FIG. 1. One-phonon correction to the fermion propagator, A phonon with momentum 2pp + k 
(dashed line) connects electron states with momenta pp +p and — pp + p — k, respectively energies 
pv and — (p — k)v. 



B. Continuum formulation 

In the vicinity of the two Fermi points, the equations of motion for the fermion field Ck(t) 
in the presence of longitudinal phonons is 

f dp 

(id/dt - vq)c PF+q {t) = gL~ 1/2 / — 5x 2 p F+q - P c_ PF+p (t), 

J Z7T 

/dp 
— Sx_2 PP + q -p c PF+p (t). (2) 

Here Sx is the atomic displacement along the chain direction, g is the electron-phonon 
coupling and L is the chain length. The electron energy spectrum near ±Pf has been 
linearized, e± PF+q ~ ±vq. 

It is convenient to combine right and left-moving fermion fields into a column 
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** = (i*)=(y +9 . )• (3) 



c -PF+q 



Phonons can be described as a complex gap field A(t,x) defined in terms of the Fourier 
transform of the displacement 

Afc(t) = g 5x 2pF+ k(t), A* k (t) = g* 5x* 2pF+k (t) = g*x^ 2pF _ k (t). (4) 

The last equality is a statement that atomic diaplacements are real. In contrast, A* k ^ A_ k , 
whereby A*(t,x) ^ A(t,x), i.e., the gap field A(t, x) is genuinely complex, except when the 
CDW is commensurate with the lattice, 2pp = it /a or 2ir/a. 
In the new notation, Eqs. (§) can be written as 

i(d/dt + v d/dx) ipR(t, x) = A(t, x) ipL(t, x), 

i{d/dt-vd/dx)ip L (t,x) = A*(t,x)ip R (t,x). (5) 

In what follows, units in which % = v = 1 will often be employed to simplify the notation. 

The fermion propagator for the ground state |0) can be defined as a 2x2 matrix G with 
matrix elements 

G aa/ (t-t',x-x') = -2(0|T[^(t,x)^,(^x')]|0>, (6) 

a = 1 for right and —1 for left fermions. Thermal Green's functions can be defined in a 
similar way. The propagator matrix satisfies the equation 



G(t-t',x-x') = S(t-t')5(x-x'), (7) 



where A(t, x) is the off-diagonal matrix A(t, x)a + +A*(t, x)<j_ and <jj are the Pauli matrices. 

The free [A(t, x) = 0] propagator G^ is diagonal in the basis of left and right-moving 
fermions, where 03 = o = ±1: 

G(°> 1/271 



v t — ox — iO sign(t) 
We will use extensively its Fourier transforms, 

G^, p) = 1 , (9) 

uj — op + i\) sign(c<j) 

G£}(u}, x) = -i sign(u;) 6(oux) e ia " x , (10) 

where 8(x) is the unit step-function. Eq. (|ToD indicates that fermions can only propagate in 
a single direction. Unless specified otherwise, it will be assumed throughout the paper that 
u > 0. 

The gap field A(t,x) is considered to be static, A(a;). In thermal field theory, this 
corresponds to a classical approximation, in which the typical frequency of a boson is much 
less than the temperature (and the occupation number of that mode greatly exceeds 1). 
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As long as this does not lead to an ultraviolet catastrophe, it appears to be a reasonable 
approximation. 

All we need now to determine the properties of the fermions are the correlation functions 
for the gap field. In the symmetric phase (above Tp), 

(A(x)> = (A*(x)) = 0, (A{x)A{x')) = (A*(x)A*(x')> = 0. (11) 

The two-point correlation function and its Fourier transform are 

D(x - x') = (A(x)AV)> = 5 2 e-\ x - x 'UZ (12) 

m = f^Tr (13) 



C. Fermion spectrum to order S 2 

The free (A = 0) fermion density of states J\f(°'(uj) can be read off directly from the 
propagator G^(uj, x): 

AT (0) M = ~ Tr ImG (0) (cu,0) = ~ £ lmG^(u,x) (14) 



7T 7T . i 

<7 = ±1 



While the value of the Green's function (fl0|) is not defined at x — 0, we can either take the 
limit x — > + or integrate over momenta 

-7r- 1 ImGW(^,p) = 5(^-ap). (15) 

Either way, the free density of states (per spin) is 

AT (0) H = 1/tt = 1/ttv, (16) 

as one expects in one dimension. 

The fermion Green's function in the presence of a gap field can be obtained by starting 
with the free propagator and iterating Eq. (|7]). This procedure gives an expansion of G in 
powers of the gap field, 

(G) = G (0) + (G (2) ) + (G (4) > + ... (17) 

The brackets signify averaging over configurations of the phonon field. The lowest-order 
correction (G^) is a diagonal matrix. E.g., for right-moving fermions, 

G { $ R {x\x) = J d(d(> G$ R (x' - C) A(C) 4°i(C - CO A*(0 GgkC - x). (18) 

Averaging over a phonon ensemble with a mean fluctuation (|I2" ) brings out (A(£) A*(£')) = 
5 2 e |C-m. Then 

(cg^, x)> = 5 2 / rfC < e" 1 ^ Gf R {x> - Q G^(C - G§»(C' - (19) 
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which is translationally invariant. 

Surely this correction can be computed in an easier way, by working directly with Fourier 
transformed quantities, Fig. |I] and Eq. ([[]). Coordinate representation, neverthelss, is also 
useful. After all, the local density of states is given by — 7r _1 Tr lmG(u; x, x). As a bonus, 
we will see where and why Sadovskii's exact solution actually works — see Sec. |111 (J| . 



"time" 




space 



FIG. 2. Second-order correction to the fermion propagator Solid lines: free fermion 

propagator —ie lujln . Dashed line: two-point phonon correlation 5 2 e~^~ ^'1/^ = 8 2 e~ l1 ^. The "time" 
direction is added to split apart fermion lines. 



Thanks to the presence of a step-function in the free propagators (10), it is more conve- 
nient to integrate over path lengths h than over intermediate coordinates Ci, C2 (Fig. ||). 
When u > 0, the free fermion propagator is 



-le 



iujln 



< l n < 00. 



(20) 



The lengths of fermion legs are not completely independent as the total displacement x' — x 
is fixed. This constraint is implemented by inserting 



5(x — x' — l + li — l 2 ) 



dP ip(x-x'-l +h-l 2 ) 

2tt 



(21) 



in the integrand. 

In particular, when we are interested in the local density of states, it makes sense to 
evaluate 



dp 



d 2 



1 



1 



dp 1 

2n uj — p + iO uj + p + i^~ x uj — p + iO 



(22) 



The integrand on the last line is the familiar second-order self-energy ([!]) with the external 
legs reattached. After the integration, we find the density of states to order 5 2 : 



W 25 2 
7iv{ (2u + iv£- 1 ) 2 



(23) 



(the factor of 2 comes from adding the contribution of left-moving fermions). The density 
of states is reduced in the range < v(,~ 1 /2, which can be called a pseudogap. This 
approximation is valid only when the fluctuations are fast enough, 5 v^ 1 /2. 



5 



D. Beyond S 2 

Correction to the fermion propagator of order S 2n reads 

(G$8&, *)> = J d( n dC n ■ ■ ■ d(x d([ G&tf - Cn) • • • Gg(Ci - CO - x) 

xD{( n ,C n ,---,(l,Ci), (24) 

where the 2n-point correlation functions is 

D(Cn, C • • • , G, CO = (A(Cn) A*(0 . . . A(C0 A*(C0) (25) 

In principle, the 2n-point correlation function must be determined in a microscopic theory. 
In most cases, however, evaluation of higher-order phonon correlation functions is a rather 
difficult job. Alternatively, one can try to see what comes out of (|24]) given a certain statistics 
of the phonon field ( ^5| ) . 

A rather trivial example would be that of the mean-field approximation, in which the 
displacement amplitude is uniform throughout the chain, i.e., A(C) =5 with certainty for 
any C- In this case, (p5| ) reduces to 

D(C„,C...,Ci,Ci)=^ 2n - (26) 

Then, the Fourier transform of the right-moving propagator is 

oo oo s2n 

(Grr(p)) = E(4 2 ?(P)) = E 7 ,! + u T w (27) 

us + pv 



oo 2 — p 2 v 2 — 5 2 
The spectral function contains two narrow peaks, 

A(u,p) = -n^lmG RR (u + i0,p) 

e p +pv sf „ ^ i e p -pu 



(28) 



^-eA + ^-^^ + g, (29) 



where 



e P = + <5 2 (30) 



(recall that p is the distance to the Fermi momentum dbp_p). The density of states vanishes 
for |a?| < S and exhibits a pile-up near u = ±5: 

M M = M -F==- (31) 

To describe a state without a long-range order, one can consider a phonon statistics with 
a fixed gap amplitude 5 and a fluctuating phase (the nonlinear 0(2) cr model). The model 
is characterized by a single parameter, a temperature-independent phase stifness a. The 
two-point phonon correlation function has been calculated, e.g., in 0: 
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D(x, x') = 6 2 e~ lx - x '^, £ = a/T. (32) 
Higher-order correlations can be computed in a similar way: 

D(x n ,x' n , . . .,xi,x[) = (A(x n ) A*(x' n ) . . . A(xi) A(x[)) 

"£ 2£ j- (33) 

As the temperature approaches zero, £ — ► oo and one recovers the statistics of the mean- 
field theory (P6f). Accordingly, the fermion energy spectrum in the limit of a long correlation 
length approaches the BCS form (|31|), something one rightfully expects. 



II. MODEL OF SADOVSKII 

A phonon system with different statistics was considered in the 1970's by M. V. Sadovskii. 
Instead of a fixed gap amplitude and Gaussian phase fluctuations (as in Sec. |1D| ), his model 
is concerned with independent Gaussian fluctuations of real and imaginary parts of A(x). 
In other words, both phase and amplitude of the gap are allowed to fluctuate. This feature 
leads to a very different fermion spectrum in the limit of slow fluctuations. 



A. Phonons with Gaussian statistics 

Statistical properties of a Gaussian random variable A(x) are completely determined 
once the mean value and two-point correlations are specified: 

(A(x)) = (A*(x)) = 0, (34) 
(A(x)A(x')} = (A*(x)A*(x')) = 0, (35) 
(A(x)A*(x')) = D(x - x') = 5 2 e-\ x ~ x '^. (36) 



All higher-order correlations fl25|) are then given by Wick's theorem, 



D(x n , x' n ,...,xi, x[) = D(x n - x' n ) ... D(x 2 - x' 2 ) D(x x - x[) 

+ permutations of primed coordinates. (37) 

The right-hand side includes n\ terms, e.g., 

D(x 2 , x' 2 , xi,x[) = 8 A e-\ X2 - x '^e-\ xl - x '^ + ^e-^-^'^e -1 * 1- ^ 1 ^. (38) 



B. Solution for £ —* oo 

As first noted by Sadovskii 0, the determination of the electron energy spectrum sim- 
plifies in the limit of long-range (slow) fluctuations of the order parameter, £ ^> 5/v . The 
electron energy spectrum in this limit is strikingly different from the spectrum with a sharp 
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gap discussed in Sec. |[D| . Instead, one finds a broadly smeared gap, or pseudogap, which is 
caused by fluctuations of the gap amplitude, absent in the previous model. 

Taking £ — > oo reduces fl36|) to a constant and thus ( |37| ) becomes coordinate-independent 
as well: 

D(x n , x' n , . . . , xi, x[) = n\ 5 2n . (39) 

By comparing this result to what we had for the state with a fixed A (p6|), we can write the 
following expression for the fermion propagator 

oo oo s2n 
(Grr(p)) = E(G'S ) (P)) = E nl 7 \~~l77 v-. (40) 

The only difference from fl27p is the factor n\, which makes the sum divergent for any fre- 
quency and momentum. 

This difficulty can be circumvented |2[ if we recognize that the divergent sum is an 
asymptotic expansion of the Stiltjes integral: 

1 hJj n=0 

This is precisely our series. The left-hand side is perfectly finite for any x away from the 
positive real axis. If x approaches the real axis from the complex plane, x ± iO, the integral 
(|41|) has a non-zero imaginary part, not reproducible by a sum of positive numbers on the 
right-hand side, hence a divergence. 

Rather than trying to resum a divergent series, it is more useful to remove the divergence 
all together. For that purpose, we will go back to the original assumption about the Gaussian 
statistics, which is the source of n\. In the limit £ — > oo, instead of a random field A (a;), 
we have a single random variable A describing the value of the gap field everywhere on the 
chain. Its Gaussian character, postulated above, is realized by considering an ensemble of 
chains, each with a different but fixed A, with the distribution ("density of chains") 

p(A) =e-^l &2 /-k5\ (42) 

which gives, as required, 

(A"A* n ) = j | A| 2n p(A) d 2 A = n\ 5 2n . (43) 

On every single chain, there is a perfect Peierls gap of size |A|, which, however, varies from 
chain to chain. 

To obtain, e.g., the density of states in such an ensemble, one can average the result for 
a single gap fl3lD over the distribution of gaps ([42]) , which gives a smeared-out gap @. The 



density of states vanishes as uo 2 at low frequencies. The fermion spectral function can be 
obtained in a similar way, by integrating the BCS spectral weight with two 5-functions ( p9[) 
over the gap distribution (fH). As a result of a varying gap amplitude | A|, one finds 0] peaks 
that are significantly broad, especially near the Fermi points (p = 0), where the spectral line 
shape is 

A(u,0) = \Lu\5- 2 e-" 2 ' s \ (44) 

A large linewidth, of order 6, reflects not a scattering rate, but rather an inhomogeneous 
broadening due to a varying gap size. 
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C. What do Sadovskii's chains look like? 



A typical Sadovskii's chain is shown in Fig. [| The gap field, or the complex amplitude 
of atomic displacements, remains approximately constant over distances smaller than the 
correlation length £. Even though this chain looks very rough, its two-point correlation 
function D(x,x') is quite smooth. In fact, 



D(X,X') = J2 e -|x-x'|/e 



(45) 



exactly for this particular chain. 

But wait. How is a correlation function defined for a single chain? And why is the 
correlation function translation invariant, 



D(x + (,x' + = D(x,x'), 
whereas the chain is not? Answer: 

D(x,x') ^ j [ L d(A(x + 0A*(x' + (). 
L Jo 

This is obviously translation invariant for periodic boundary conditions. 



(46) 



(47) 




FIG. 3. Half of a typical Sadovskii chain. The correlation length (the two-head arrow) is £ = 1/8, 
the mean value of the gap is 5 = 1. Solid line: ReA(x), dashed line: ImA(x). The actual atomic 
displacement at a point x is given by Re A(x)e 2lkpX . 

To see how a chain with the right correlation function can be constructed, rewrite ( fT7| ) 
in terms of Fourier components 



D(x,x') = L- 1 f L d( L~ l l 2 Y,^ke lkix+0 L-W'EXve-'W- 
Jo k k> 

= L- 1 \A k \ 2 e lk{x ~ x ' ] ^ \y \A k \ 2 e ik{x ~ x ' ] 



(4* 
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in the limit L — > oo. Choosing now 



with an arbitrary phase 9k immediately yields (|45|) . This is approximately how the chain in 
Fig. [3] has been simulated. 

Furthermore, it can now be seen that an ensemble of such chains exhibits the required 
Gaussian statistics (]5T|) . For instance, the four-point correlation function for a single chain 
is 

D(x 2 , x' 2 , x 1 ,x' 1 ) ^ \ [ L d( A(x 2 + C) A(x' 2 + C)A(x! + C) A*(x[ + () 
L Jo 

= L- 2 &k 2 e ik2X2 A* k ! 2 e- lk! ^ A kl e lklXl A* K e~ iK< , (50) 
{k} 

where momenta satisfy the constraint k\ + k 2 = k[ + k' 2 . 

The crucial step is to average ( ^0[) over the arbitrary phases 9k, which enter this expression 
in the form of the factor 



exp 



i{6 kl +9 k2 -9 k , -9u) . (51) 



Phases 9k are independent random variables uniformly distributed in the interval < 9 k < 
2tt. Averaging over them makes (^Tj) vanish, unless the phase factors cancel one another 
pairwise: 

ki = k[, k 2 = k' 2 , or k\ = k' 2 , k 2 = k[. (52) 
The so averaged four-point correlation function reads 

(D(x 2 ,x' 2 ,x 1 ,x' 1 )) = L-^l^k^e^ 2 '^ L-^IAfcjV' 11 ^ 1 -^ 

k 2 ki 

+ L- 1 Y,\&k 2 \ 2 e ik2{x2 - x 'i ) L-^IAfcjV^ 1 ^ 1 ^) (53) 
fc 2 fci 

_ ^^\i^ik(x 2 -x' 2 +x 1 -x' 1 ) 
k 

The third line in (|53|) is needed to adjust for the overcounting of the terms with k\ = k[ = 
k 2 = k! 2 in (f)2|). Comparing these to (fP|) reveals (almost) Gaussian statistics: 

(D(x 2 , x 2 , xi, x\)) = D(x 2 , x' 2 ) D(x 1 , x\) + D(x 2 , x[) D(x 1 , x' 2 ) (54) 

-IT 1 l—\[\ k \ A e ik ^- x '^^-0_ 
J 2n 

In the limit L — > oo, the extra term vanishes as £/L or faster. This procedure can evidently 
be extended to arbitrarily high orders. 
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FIG. 4. Four-point correlation functions (D(x,x,0,0)) (crosses) and {D(x, 0, x, 0)) (diamonds) 
averaged over 32 different chains. Solid lines are the corresponding Gaussian curves. £ = L/64, 
5 = 1. 

We thus have found an efficient way to simulate the ensemble considered by Sadovskii. 
Namely, by using the known Fourier amplitudes (|Hf) , we can generate a sufficiently large 
number of chains with different choices of random phases 9k- That this number does not have 
to be too large is demonstrated in Fig. [|. Average four-point correlations in an ensemble of 
just 32 chains agree quite well with their Gaussian expectation values. Note a systematic 
downward shift for (D(x, x, 0, 0)), which is expected to equal £/L. 

As a matter of fact, it is not necessary to consider an ensemble of chains. Instead, one 
can regard a single long enough chain as an ensemble of its segments (of length L) . Doing so 
resolves the apparent paradox of strongly broadened electron states in the limit £ — > oo. As 
long as one studies properties of electrons within a single segment, this broadening will exist 
until £ exceeds the segment length L. In the limit £ 3> L, electrons will see a well-defined gap 
A (within this segment) and their spectral function will exhibit two sharp Bogoliubov peaks 
at energies ±yp 2 t> 2 + |A| 2 . However, if we now pack our instruments and go to another 
segment of the chain, far enough away, we could find there a different value of the gap. 
Averaging over an infinitely long chain, on the other hand, will always give the Gaussian 
broadening as £ < L = oo always in that case. 

Remark. The possibility of gap-amplitude fluctuations is largely ignored in the literature. 
Common-sense wisdom suggests that variations of the energy gap cost too much energy 
and therefore only the phase of the order parameter A is expected to fluctuate at low 
temperatures. There are two counterarguments here. First, a qualitative one, that gap- 
size fluctuations further increase the entropy and thus reduce the free energy. Another, 
quantitative argument appeals to the well-known solution || for the fermion spectrum in the 
Luttinger model with attraction, truly an exactly solvable model with a pseudogap. Cooper 
pairs, the dominating fluctuations in this model H, produce a broadly smeared energy gap 
in the fermion spectrum, even at zero temperature. In the limit of a large correlation length, 
phase fluctuations alone cannot account for a strong smearing. Thus, amplitude fluctuations 
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are inevitably present in this model and may be more commonplace than usually thought. 



III. SADOVSKII'S SOLUTION 

We now return to the theoretical analysis of the problem, this time for a finite correlation 
length £. 

A. Sadovskii's conjecture 

As previously noted, evaluation of Feynman diagrams is more convenient in momentum 
space, where the two-point phonon correlation is a Lorentzian (|T3D, 



D (k) = S 2 J^- 2 (55) 
As we have seen before, the order-5 2 correction to the electron Green's function is (Fig. [[]) 
G W M = l —— . (56) 

This line is rather transparent: integration over momentum k transfered to the phonon 
simply shifts the imaginary in the denominator of the intermediate electron propagator from 
+i0 to -H£ -1 . 

(a) (b) /" \ 

— • — • — ■< m * — 

FIG. 5. Two-phonon contributions to the fermion propagator, (a) is generated by the first-order 
self-energy, while (b) contains two-phonon self-energy. 

Such a simple form of the second-order correction, 

sgi^^cxG^^ + ^r 1 ^), (57) 

has prompted Sadovskii to conjecture that contributions of higher-order graphs to G a<T (p) 
are given by the following simple rules: 

A phonon line contributes 5 2 , (58) 

An electron line contributes — — , (59) 

u±p + iv^- 1 y ' 

where v is the number of phonon lines above a given electron line. The sign in front of p 
alternates as the fermion propagates left and right. 

For instance, according to this rule, corrections of order 5 A to the electron Green's function 
(Fig. |5|) should read 
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co> — p + iO I u — p + 2i£ 1 V to + p + i£ 1 



(60) 
(61) 



Basing on this Ansatz, Sadovskii was able to derive and solve a recursion relation for the 
self-energy of order 5 2n [|3| following a method due to Elyutin |J . The exact Green's function 
was then obtained in a continued fraction representation. This remarkable derivation is 
getting quite popular these days (|,[J. 



B. Failure in order <5 4 

Unfortunately, Ansatz (^9|) works only for a limited class of diagrams [e.g., Fig. |5](a)] and 
is simply incorrect for others [Fig. |5|(b)]. The problem, quite mundane, is in sloppy handling 
of the imaginary part — equal to +i0 sign(u;) or iuo n depending on the formalism — in the 
denominator of G^(uj,p). 

Recall that, to order 5 2 , we integrated 

dk 1 2^ 

2itu + p-k + iQ P + £" 2 ' 1 ' 

which has two poles above the real k axis and only one pole below. If we complete the 
integration contour in the lower half of the plane, only one pole is inside and the resulting 
expression is simple. 

In evaluating graph (b) in Fig. |5], the integral over momentum q of the external phonon 
line reads 

dq 1 1 m 

2tt [u + (p-q)+ iO} 2 uu — (p — q-k)+i0 q 2 + £~ 2 ' 1 ' 

which has two poles on either side of the real q axis, so that, whichever way the contour is 
completed at infinity, the result contains two terms, rather than one. Integrating over k first 
does not help either: 

* I ! ?H- (64 ) 

2tt [uj + (p - q) + iO} 2 uj - {p - q) + q 2 + £~ 2 V ; 

is plagued by the same problem. The integral over q, completed above the real axis, yields 
the result conjectured by Sadovskii ( |6"TD plus a non-zero contribution from the pole at q — 
uj + p + iO. 

The situation does not change when one uses thermal Green's functions, in which case 
u + iO sign(cu) is replaced with iuj n and the same problem arises. In higher orders, expressions 
for G^(p) becomes progressively more complicated by the presence of diagrams with a 
phonon line running over many electron propagators. 
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C. Which problem did Sadovskii solve, exactly? 



(a) 




"time" 



^2 Ci 



(b) 




C 2 Ci 



FIG. 6. Correction to the fermion propagator Solid lines: free fermion propagator 

— ie lwln . Dashed lines: two-point phonon correlation <5 2 e _ '^ i- ^ . Vertical dimension is added for 
clarity. 



The fact that the trouble is caused by infinitesimal imaginary numbers in fermion propa- 
gators may create an illusion that the problem can be somehow fixed. It is more instructive 
to look at it in configuration space. We will now see exactly which problem Sadovskii solved. 

Using conventions of Sec. [TC], we write out the expression for the first of the two diagrams 
for (Gi&(p)>, Fig. |(a): 



poo 



x j4 g-lfia-^l/C e -|0-Cil/5 



(65) 



[cf. Eq.fl22|)]. As IC2 — C2I = ^3 an d ICi — Cx I = hi the integrals over lengths {/„} can be 
immediately carried out and one obtains (|5T)|). 

The other diagram, Fig. §(b), differs by a permutation of ([ and so that only the 
second line of (BSJ) changes and now reads 



x 5 4 e" IC2 "^ l/5 e ~ ICl ~ ? 2 |/f . 



(66) 



While |Ci — C2I = hi the other distance, \( 2 — Cil; cannot be simply expressed as a sum of 
some path lengths, which is what causes the problem. Note, however, that, had we replaced 
the physical distance IC2 — Cil with the sum of path lengths h+h + h, the previous expression 
would have read 



x S 4 e~ 1 ^ e~ 21 ^ e^, 



(67) 



which could be easily integrated over lengths yielding Eq. Q5"T|), precisely what Sadovskii 
wanted. 
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Once the physical distance between two points \Q — £j| in the phonon correlation func- 
tion D(d,(j) has been replaced with the length of the fermion path between these points, 
Sadovskii's conjecture fl56|) is valid in all orders of perturbation theory. Indeed, define v m to 
be the number of phonon lines above the fermion leg l m , which can be done unamiguously 
by straightening out the fermion trajectory (i.e., by using Fig. [5] instead of Fig. ^]). The 
contribution of a given diagram to G 2 ^ R {p) will then be a product of independent factors 



^ n />00 



m=0 

2n i 

= I> J1 »-Hh^ ' (68) 



precisely as required by 

It is thus clear that the original Ansatz of Sadovskii solves a rather unphysical problem, 
in which phonon correlations (A(x)A*(x')) depend not on the geometrical distance \x — x'\, 
but rather on the length of the path the fermion traveled between points x and x' . This 
point is further illustrated using a two-dimensional example in Sec. [TV| 

IV. EXTENSION TO HIGHER DIMENSIONS? 

It has already been mentioned that the calculation of Sadovskii is tailored to one spatial 
dimension. This limitation stems from the fact that the order-5 2 correction to the fermion 
self-energy in the presence of classical fluctuations of an order parameter, 

1 (69) 



(2n) d u + (p-k)-v + iO fc 2 + £~ 2 ' 

has a simple form in d = 1 dimension only. Not having a simple, "self-replicating" form for 
the lowest-order correction possibly indicates that there is little hope of finding a general 
recipe for higher orders. 



A. Model of "hot spots" in the cuprates 

Recently, however, J. Schmalian, D. Pines and B. Stojkovic applied the ideas of 
Sadovskii to a two-dimensional system, high- Tp cuprate superconductors, to investigate the 
nearly antiferromagnetic Fermi liquid (NAFL) flQ| . This development further illustrates in 



what context the solution of Sadovskii is applicable. It turns out that the dimensionality 
of the system is not important. A really necessary ingredient is the peculiar form of order- 
parameter correlations, which should decay exponentially with the "distance" measured 
along the fermion path. (There is also a technical, but very important, requirement that the 
spectrum of free electrons be flat, i.e., e PF+p = p-v, where v is a constant vector.) 
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71 r' 

(a) (b) 

FIG. 7. (a) A sketch of the Fermi surface (solid line) in the cuprates. vi and v 2 are Fermi 
velocities at two hot spots (filled circles) connected by the antiferromagnetic wave vector Q = (it, tt). 
The dashed line is the locus of states most strongly affected by the AFM scattering, e p+ Q = e p . 
(b) A fermion initially in the vicinity of the hot spot pi travels in a zigzag manner in real space 
switching between non-collinear velocities Vi and V2 as it is scattered by spin fluctuations. 

In the NAFL approach, electrons, considered to be ideal fermions, interact with antifer- 
romagnetic (AFM) spin fluctuations, whose static susceptibility is peaked near wavenumber 
Q = (tt, 7r) in reciprocal lattice units: 

x(Q + q)-7^- (70) 

The strongest effect of AFM fluctuations on the fermion energy spectrum is expected when 
scattering by wave vector Q connects states of the same energy, e p+ Q = e p . Such points in the 
Brillouin zone form a line shown in Fig. ^(a) for a tight-binding fermion energy spectrum 
(nearest and next-nearest neighbor hopping). Places where this line intersects the Fermi 
surface have been termed "hot spots". Low-energy fermionic excitations in these spots are 
presumably fried by spin fluctuations and are short-lived, hence the name. This must be 
true, at least to some extent, as photoemission shows extremely broad peaks (hundreds of 
meV) in the electron spectral weight A(uj, p) at these momenta |pT |. 



Consider the lowest-order fermion self-energy from one-magnon exchange (the same dia- 
gram as in Fig. [I]). After linearizing the free fermion spectrum near the hot spots 



L Pn+P 

the self-energy for a fermion near hot spot pi reads 



~v„-p, (71) 



^■p. + pWt^t^U t r wn - (72) 

J {2-k) 2 1 + <t4 2 uj + v 2 -(p - q) + lO 

Here Q + q is the momentum transfered to the magnon. Note that the intermediate electron 
is near the other hot spot p 2 . This has pecisely the form of Eq. ( |69|) and one cannot get a 
simple expression out of it, to say nothing of higher-order corrections. 
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Schmalian et al. noted that Fermi velocities at conjugated hot spots (e.g., vi and V2) 
are almost perpendicular to each other. If then one replaces the susceptibility (|70|) with a 
product 



x(Q + q)-5 2 ^-I-^-^, (73) 



2£ _1 2£" x 

where q n is the component of q along v n , a very simple sef-energy results: 



E(u;, Pl + p) « — 6 \. cl = (u + ivC\ P2 + P), (74) 

u + v 2 p + 1 

where v = |vi| = |v 2 |. Moreover, higher-order diagrams can be evaluated in a similar manner 
yielding simple expressions in the form conjectured by Sadovskii (|56]) . Lo and behold, the 
problem becomes tractable to arbitrary order and the electron Green's function can be 
obtained in the continued fraction representation H, as discussed by Sadovskii. 



B. What makes it solvable 

One should not be surprised that the trick with factorization (|73|) makes the problem 
solvable. The factorization amounts to taking spin-spin correlations in real space in the form 

X (r - r') oc (s + (r) S -(r')> oc e-^'^e^'^, (75) 

where x and y are Cartesian components of the electron in the plane (along the directions 
of vi and v 2 , i.e., approximately along the crystal axes). This is precisely the same as to 
say that order parameter correlations decay with the "distance" measured along the fermion 
path, Fig. ^(b), as I noted in the beginning of this Section. That said, it is not even necessary 
to require that vi and v 2 be orthogonal. 

Technically, the similarity with the one-dimensional problem of Sadovskii arises because 
the electron energy spectrum has been linearized, whatever the actual number of dimensions 
is. Indeed, according to (|7T|) , the energy as a function of momentum varies only in the 
direction of v n . The fermion spectrum is exactly dispersionless in all directions perpendicular 
to v n . This means, literally, that a free fermion propagates along a straight line, Fig. 0(b). 
In this approximation, d 2 e/dpidpj = 0, a localized wave packet does not disperse as it 
propagates with velocity Vi = de/dpi. It would not be a stretch to say that this is essentially 
a one- dimensional problem. 



C. No pseudogap in the DOS 



Despite great similarities, there is one important aspect in which this two-dimensional 
problem differs from the purely one- dimensional case of Sadovskii. In plain English, a fermion 
never returns to a starting point: it zigzags away, Fig. 0(b). In contrast, Vi = — v 2 in one 
dimension and a fermion does return to the starting point "once in a while". An impor- 
tant consequence of this innocuous observation is that, with vj 7^ — v 2 , the local fermion 
propagator is unaffected by fluctuations to all orders, 
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G(u,x, x) = G^°\uj,x,x), 



(76) 



because (u, x, x) = for any order n > 0. Therefore, the local density of electron states 
is exactly the same as for free, noninteracting fermions: 

Af{u) = -Tr^ImG^^x) = -T:- 1 lmG^\uj,x,x) =AT {0) {uj). (77) 

In particular, this means that a local probe, such as tunneling microscopy |T3 or NMR [O 



should not observe any pseudogap behavior in the "hot spots" scenario! 

The use of the coordinate representation makes a proof of this statement almost trivial: 
the step function in the free electron propagator (|T0| ) makes all corrections to the free prop- 
agator vanish for fermion trajectories with points outside the cone formed by the vectors Vi 
and v 2 . Thereby trajectories returning to the starting point consist of a single point, have 
zero integration measure (in the case of at least one intermediate point) and therefore do 
not contribute to the propagator. 

Of course, this can be seen in momentum space as well. For simplicity, take the velocities 
Vi and V2 to be orthogonal to each other and choose a pair of coordinate axes along them. 
A generic correction of order S 2n to the propagator of a fermion near Hot Spot 1 contains 
two momentum-dependent factors: 



where v\ a ,V2b > (strict inequality!). Its contribution to the local DOS is obtained by 
integrating over the momentum components p\ and p 2 and taking the imaginary part. This 
expression is an analytical function of p\ below the real axis. If n > 0, the integration 
contour can be completed at infinity in the lower half of the complex p\ plane (the integrand 
vanishes there fast enough). Since no singularities are encircled, the integral vanishes for 



any n > 0. The exceptional case n = (free propagator) has been dealt with in Sec. |IC . 

It is worth stressing that Sadovskii's solution should be considered as a long-wavelength 
approximation only (as it is based on a linearized electron spectrum). In practical terms, 
one should not attempt to draw conclusions about the detailed band structure basing on a 
solution of this type. Things like a Brillouin zone or a van Hove singularity simply do not 
belong in this theory. While integrating the spectral weight over a Brillouin zone may show 
a slight reduction in the density of states near the Fermi level ||14|| , such an extrapolation of 
an effective field theory to real-life details is not warranted. The only conclusion that can 
be drawn safely is that, as the lattice spacing is taken to zero, any trace of the pseudogap 
disappears. Thus, the pseudogap is not natural in this model. 
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(0,pi) hot spot 



(pi/2,pi/2) 
wavenumber 



hot spot (pi,0) 



FIG. 8. Splitting of the free electron band (dashed line) into upper and lower bands (solid lines) 
in the presence of long-range AFM order with wavevector Q = (ir/a, it /a). When the band splitting 
25 is smaller than the bandwidth, there is no gap in the density of states. 

While it may appear paradoxic that the DOS is unaffected, it is, in fact, a direct con- 
sequence of the assumptions that made the calculation of Schmalian et al. possible. It is 



also directly related to an observation by Randeria |15| that the NAFL pseudogap is not 
tied to the Fermi surface (the dotted line and the solid line in Fig. ^|(a), respectively). In 
the antiferromagnetic scenario for the pseudogap, the spectral weight of the fermion states 
on the dotted line is moved from e p to higher and lower energies in the range e p ± 5. Since, 
however, the energy e p varies along the dotted line (by the amount equal to 4i 2 , where t 2 is 
the next-nearest neighbor hopping amplitude), the pseudogap will be completely washed out 
if 4t 2 exceeds 6. Linearization of spectrum ([H]) is equivalent to assuming 5 <C 4t 2 (no local 



pseudogap). As Monthoux and Pines suggested, At 2 = 0.45 eV JTIJ, so that any pseudogap 
of a lesser width will be washed out in the density of states. This situation is illustrated 
schematically in Fig. [8]. 

In contrast, there is no washing out of a pseudogap created by Cooper pair fluctuations. 
In that case, fermion states coupled by emission or absorption of a Cooper pair are electrons 
and holes of equal momenta and spin. Therefore, their velocities are equal and opposite, 

vji= rfh !zE ) = _* E = 

dp dp 

as long as time reversal is a good symmetry of the system (e_ p = e p ). In the problem with 
a linearized dispersion, a fermion moves along a straight line back and forth alternating 
between an electron and a hole. There are non-zero corrections to the local propagator 
G(u,x, x) in all orders, which means that scattering by Cooper pair fluctuations does affect 
the local DOS. Put simply, a pseudogap created by pairing fluctuations is tied to a Fermi 
surface. 
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V. SUMMARY 



In this paper, it has been demonstrated that the issue of phonon statistics is quite impor- 
tant for the properties of electrons in the pseudogap regime above the ordering temperature. 
Knowledge of the two-point correlation function (A(sc) A* (a;')) allows one to compute the 
electron Green's function or self-energy to the second order in the gap size 5 only. When 
the correlation length of the fluctuations increases beyond the point £ > v/8, higher-order 
phonon contributions become important, which is why multi-phonon correlation functions 
are needed. It has been shown explicitly that different choices of phonon statistics lead to 
widely different results for the fermion spectrum in the particularly interesting limit of slow 
fluctuations, £ ^> v/8. 

A model of phonons with Gaussian staistics 0j has been revisited and thoroughly dis- 
cussed, both in momentum and coordinate domains. It has been shown that its "exact" 
solution for a finite correlation length |§ contains an error and, in fact, solves another, 
rather unphysical problem. 

The physical reason why the gap in the density of states remains smeared even for very 
long correlation lengths in the model with Gaussian phonons resides with the fluctuations 
of the gap amplitude inherent in the model. This smearing should not be interpreted as a 
presence of a large (of order 5) scattering rate. Rather, it should be regarded as an inhomo- 
geneous broadening of energy levels, which, being a reversible process, can be distinguished 
from relaxational broadening. To do so, one may attempt to study the fermion lineshape 
using time-resolved spectroscopy rather than frequency-domain methods (cf. NMR). 

Finally, I have discussed a few aspects of the newly proposed scenario for the behavior 
of electrons at "hot spots" in cuprate superconductors ||. In particular, it appears that the 
antiferromagnetic fluctuations alone cannot explain the presence of a strong pseudogap seen 
by local probes of the density of states, such as tunneling spectroscopy and NMR. As has 
been noted before ||15|| , pairing fluctuations seem to be a necessary ingredient to explain the 
pseudogap at low frequencies. 
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